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We apply the derivative expansion approach to the Casimir effect for a real scalar field in d 
spatial dimensions, to calculate the next to leading order term in that expansion, namely, the first 
correction to the proximity force approximation. The field satisfies either Dirichlet or Neumann 
boundary conditions on two static mirrors, one of them flat and the other gently curved. We show 
that, for Dirichlet boundary conditions, the next to leading order term in the Casimir energy is of 
quadratic order in derivatives, regardless of the number of dimensions. Therefore it is local, and 
i—{ determined by a single coefficient. We show that the same holds true, if d 7^ 2, for a field which 

satisfies Neumann conditions. When d = 2, the next to leading order term becomes nonlocal in 
coordinate space, a manifestation of the existence of a gaplcss excitation (which do exist also for 
1 d > 2, but produce sub-leading terms). 

We also consider a derivative expansion approach including thermal fluctuations of the scalar field. 
We show that, for Dirichlet mirrors, the next to leading order term in the free energy is also local 
for any temperature T. Besides, it interpolates between the proper limits: when T — > it tends to 
the one we had calculated for the Casimir energy in d dimensions, while for T — > 00 it corresponds 
to the one for a theory in d — 1 dimensions, because of the expected dimensional reduction at high 
temperatures. For Neumann mirrors in d = 3, we find a nonlocal next to leading order term for any 
T > 0. 

Q_l PACS numbers: 12.20.Ds, 03.70.+k, ll.10.-z 
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I. INTRODUCTION 

The determination of the Casimir force [1] for a quite general situation, namely, when the geometry of the problem 
is characterized by two rather arbitrary surfaces, is interesting and potentially useful. Those surfaces may correspond, 
for example, to the boundaries of two mirrors. Alternatively, the surfaces themselves may describe zero-width ('thin') 
mirrors, which will be the situation considered in this paper. Yet another possibility is that those surfaces may be the 
interfaces between media with different electromagnetic properties, occupying different spatial regions. In situations 
like the ones above, it may be convenient to think of the Casimir energy as a functional of the functions determining 
the surfaces. Of course, it is generally quite difficult to compute that functional for arbitrary surfaces. Rather, exact 
results are only available for highly symmetric configurations, the simplest of which being the case of two flat, infinite, 
parallel plates. Taking advantage of the simplicity of the result for this highly symmetric configuration, the proximity 
force approximation (PFA) [2] 13] provides an accurate method to calculate the Casimir energy when the surfaces are 
gently curved, almost parallel, and close to each other. Introduced by Derjaguin many years ago [5] to compute Van 
der Waals forces, this approximation consists of replacing both surfaces by a set of parallel plates. The energy is then 
calculated as the sum of the Casimir energies due to each pair of plates, each plate paired only with the nearest one 
in the other mirror. 

In a recent work [4], we have shown that the PFA can be put into the context of a derivative expansion (DE) for 
the Casimir energy, when the latter is regarded as a functional of the functions that define the shapes of the mirrors. 
Indeed, the leading order term in this expansion, which contains no derivatives, does reproduce the PFA, while the 
higher order ones account for the corrections. In that article we considered, for the sake of simplicity, a massless 
quantum scalar field satisfying Dirichlet boundary conditions on two surfaces. One of them was assumed to be flat, 
and such that if coordinates were chosen so that x 3 = 0, the other surface could be described by a single function: 

Since the form of the possible terms in the DE may be determined by dimensional analysis plus symmetry consid- 
erations, what is left is the calculation of their respective coefficients. Note that those coefficients are 'universal', in 
the sense that they are independent of the shapes of the mirrors (at least for smooth surfaces) . Therefore one can fix 
those coefficients completely, from the knowledge of their values for a particular surface, or for a family of surfaces. 
We used, to that effect, a particular family of surfaces, namely, those obtained by an expansion up to the second 
order in 77, which is the (assumed small) departure from flat parallel mirrors: ip — a + tj, r\ << a. The coefficients 
determined from this family were then used to fix the coefficients of the first two terms in the DE, which can then be 
used to calculate the Casimir energy for more general (but smooth) surfaces. 

This approach has been generalized by Bimonte et al [51 [5] in many directions. For example, to the case of two 
curved, perfectly conducting surfaces, for scalar fields satisfying Dirichlet or Neumann boundary conditions, and also 
to the electromagnetic case and imperfect boundary conditions. As a validity check, it has been shown that, whenever 
analytic results arc available for particular geometries, the corresponding DE does reproduce both the PFA and its 
next to leading order (NTLO) correction [H [5] . The DE approach has also been applied successfully to compute the 
electrostatic interaction between perfect conductors [7j. 

In [5j , we have extended our previous work 0] to the case of the electromagnetic field coupled to two thin, imperfect 
mirrors, described by means of the vacuum polarization tensors localized on the mirrors. We have also calculated 
the NTLO to the PFA static Casimir force. For the particular case of mirrors described by a single dimensionless 
quantity, we have computed the leading and NTLO corrections as a function of that quantity. We found that the 
absolute value of the NTLO correction falls down rather quickly for imperfect mirrors [S]. 

In this article, we apply the DE approach, to models where the fluctuating field is coupled to two perfect mirrors, 
L and R, at zero or finite temperature, in d + 1 spacetime dimensions. Special consideration shall be given to the 
limiting cases where the fluctuations are either thermal or quantum mechanical, i.e, infinite or zero temperature. 

The role of the fluctuating field is played by a massless real scalar field (p, with either Dirichlet or Neumann boundary 
conditions on both mirrors. This analysis is of interest because of several reasons: on the one hand, as a first step 
towards incorporating thermal effects into the DE (a fuller treatment should also include finite conductivity corrections 
along with finite temperature corrections [S]). On the other hand, as we will see, this more general analysis will lead 
us to a clearer physical picture of the validity of the DE and will also shed some light about possible extensions and 
improvements . 

We will mostly consider the first two terms in the DE; in order to fix their coefficients, we follow the procedure of 
expanding the vacuum energy up to the second order in r\ and then extracting the coefficients from the corresponding 
momentum space kernel. Thus, in what may be considered as a byproduct of our approach, we also present the 
general result for that kernel, valid for any d, both for the Dirichlet and Neumann cases. 

This kernel for the quadratic term in n (regarded now as a field) may be interpreted as a contribution to its 2-point 
one-particle irreducible function, due to a one-loop <p contribution, which fluctuates satisfying the proper boundary 
conditions. The coefficients of the DE up to the NTLO result from that kernel, from its expansion up to the second 
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order in k, the 77- field momentum, as one would do in an effective field theory approach [5]. 

When that expansion to order k 2 does exist, the NTLO terms are of second order in derivatives, and therefore 
spatially local. However, as in any quantum correction to a 2-point function, we know that non-analytic contributions 
may arise when the external momentum reaches the threshold to excite modes of the field in the loop. In the present 
case, since we expand around k — 0, those non-analyticities may only originate in the existence of massless modes. 
When present, they result in contributions which, being nonlocal in space, are similar to the ones that appear in the 
context of effective field theories, when including the effect of massless virtual particles [9]. We shall see that this 
kind of non-analyticity does indeed appear, for Neumann conditions, in the form of branch cuts. The physical reason 
being that when both mirrors impose Neumann conditions there are transverse gapless modes for the fluctuating field. 
However, except for d = 2 and zero temperature (T = 0), or d — 3 and T > 0, those non-analyticities are of higher 
order than the NTLO, when k — > 0. In this, the only 'pathological' case (d = 2 and T = or d = 3 and T > 0), the 
real time version of the kernel has, if rotated to real time, a logarithmic branch cut at zero momentum of the form 
k 2 \og(k 2 ), which overcomes the k 2 term (which is also present). 

We also study the NTLO term as a function of temperature. In the particular case of d = 3, we shall see that 
for Dirichlet boundary conditions it depends smoothly on the temperature. For Neumann boundary conditions, as 
an indirect consequence of the non-analyticity at d = 2 and zero temperature, the Neumann NTLO term is also 
non-analytic at any non zero temperature. 

This paper is organized as follows. In the next Section we introduce the system and summarize the approach 



we follow to calculate the free energy Tp. In Section III we discuss the DE for Dirichlet boundary conditions on 



thin, perfect mirrors in d + 1 dimensions, discussing the zero and high (infinite) temperature limits. We apply those 
results to evaluate the Casimir interaction energy between a sphere and a plane at very high temperature. Section |TV| 
is devoted to study the DE at zero and high temperatures limits, for a real scalar field with Neumann boundary 
conditions on the mirrors. The special cases of d = 2 with T = 0, and d = 3 with T > 0, are singled out and dealt 
with in subsection |IVB I Higher order terms in the DE are analyzed in Section|V] Finally, in Section |V*l| we summarize 
our conclusions. The Appendices contain some details of the calculations. 



II. THE SYSTEM 



We shall adopt Euclidean conventions, whereby the spacetime metric is the identity matrix, and spacetime co- 
ordinates are denoted by x^ = x^ (/i = 0, 1, . . . , d), Xq being the imaginary time and Xi, (i — l,...,d) the spatial 
Cartesian coordinates. 

Regarding the geometry of the system, we shall assume that one of the surfaces, L, is a plane, while the other, R, 
is such that it can be described by a single Monge patch: 

L) x d = R) x d = ip(x ,xi,...,x d -i) ■ (1) 

We have included for R in a more general, time-dependent boundary, in spite of the fact that we are interested 
in the static Casimir effect (SCE). We shall, indeed, at the end of the calculations, impose the condition that the 
boundaries are time-independent: ip — ij){xi, . . . ,Xd-i), but it turns out to be convenient to keep the more general 
kind of boundary condition at intermediate stages of the calculation. In this way, the treatment becomes more 
symmetric, and one may take advantage of that to simplify the calculation. Besides, although it is not our object in 
this paper, one could rotate back some of the results thus obtained for a non-static ip to real time, in order to consider 
a dynamical Casimir effect (DCE) situation. 

We follow a functional approach to calculate the free energy Tpfy), or its zero temperature limit E vac (^), the 
vacuum energy. Both are functionals of tp, that defines the shape of the R mirror (the plane mirror L is assumed to 
be fixed at Xd = 0). is also a function of the inverse temperature ft = T^ 1 (we use units such that Boltzmann 
constant ks = 1). 

In the functional approach, which we shall follow, both objects are obtained by performing a functional integration; 
indeed: 

W) = -ilog[^M]. (2) 

where Zp{tp) is the partition function; it may be obtained by integrating over field configurations that satisfy the 
corresponding boundary conditions at L and R, and are also periodic (with period (3) in the imaginary time coordinate 

xq (Matsubara formalism). Zf' denotes the partition function in the absence of the mirrors; therefore it corresponds 
to a relativistic free Bose gas. 
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-E'vac(V') is then obtained by taking the limit: 

£vacW = lim r (3 (V>) = r^ty) . (3) 

It is worth noting that in ([2]) i/j has to be time-independent for Tp(ip) to be a free energy. We do however calculate 
objects like Zp(ip) for configurations that may have a time dependence, keeping the same notation. 

To avoid an unnecessary repetition of rather similar expressions, we shall write most of the derivations within the 
context of a finite temperature system, presenting their zero temperature counterparts at the end of the calculations. 

We shall consider a real scalar field, with either Dirichlet or Neumann boundary conditions. In both cases, the 
general setup has a similar structure, but there are also some important differences. Mostly, they come from the 
different infrared behaviour of their respective Green's functions, and the impact that that behaviour has on the DE. 
Accordingly, we present them in two separate sections. 



III. DIRICHLET BOUNDARY CONDITIONS 



We start from the functional representation of Zp(ip): 

Zpty) = fvip S L (<p) 5 R (<p) e- s ^ , (4) 



where 5a(<p), A — L, R, is a functional 6 function which imposes Dirichlet boundary conditions on the respective 
mirror, while Sq is the free Euclidean action for a massless real scalar field in d+ 1 dimensions, at finite temperature: 

S o = ^J dx o Jd d x(d V >) 2 , (5) 

with periodic conditions for (p in the time-like coordinate, namely, ip(xo,x) = <p(xq +/3,x), for all x e M.( d \ 

To proceed, one should then exponentiate the 5-functionals by introducing two auxiliary fields, Xl and Xr, functions 
of X|| = (xo, xi, . . . , Xd—i) = (xq, xy), also satisfying periodic boundary conditions in the Xq coordinate. In the Dirichlet 
case, we have: 

5 L (ip) = J VX L e l S d " x " M*iiM*ii>o) 

S R ((p) = J VX R e^^n Vs( x iO x r^ x \0 ^nM^n)) ; (6) 

where g is the determinant of g a p, the induced metric on R: 

g a p(x\\) = 5 a p + d a ip(x\i)dpip(x\i) , (7) 

=> 5(x„) = l + (9V(x||)) 2 . (8) 

We have adopted the convention that indices from the beginning of the Greek alphabet run from to d — 1. 
Using the exponential representations above in Q, one derives the alternative expression: 

z p {$) = J VipVX L vx R e -so( V ) + ifd d+i xj D (x) V (x) : ^ 

where the 'Dirichlet current' Jd{x) is given by: 



Jd(x) = A l (.T||)<5(.t 3 ) + \ R (x\\) yJgixwWxs - i/j(x\\)) . (10) 

It is possible to get rid of the y/g(x\\) factor above just by redefining A^: A^(x||) — > Xr{x»)/ ^Jg{x\\). This redefinition 
induces a nontrivial Jacobian. However, this Jacobian is independent of the distance between the mirrors, therefore 
irrelevant to the calculation of their relative Casimir force; hence we discard it. 
The integral over (p, a Gaussian, yields: 

Z fi M = Zf [ VX L VX R e^H, M-i^Ci^Wi^ (n) 
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where we have introduced the objects: 

r LL {x h x\) = (a^oK-a 2 )- 1 !^) 
TThl^II,^) = (xp^x^K-af)- 1 ^) 

where we use a "bra-ket" notation to denote matrix elements of operators, and 



E 



(2tt) c 



i<4 + k2 ) 



A(ar-i/) , 



where we have introduced the Matsubara frequencies: u}„ 



27m 



n G Z. The free energy (■(/>) is then 



1 

W 



TrlogT 



(12) 
(13) 
(14) 
(15) 

(16) 
(17) 



where ip is regarded as time independent, something that one can impose at the end of the calculation. 

Tp still contains 'self-energy' contributions, i.e., contributions invariant under the rigid displacement ^(arii) — > 
Tp{x\\) + e. Since we are just interested in the Casimir force we shall neglect them altogether whenever they emerge in 
the calculations below. 



A. Derivative expansion 



The DE is implemented by following the same idea and approach introduced in [3]. The calculation is, in many 
aspects, identical to the one in [3], with the only differences in the number of dimensions and in the fact that the time 
coordinate is compact (periodic), so frequency integrations have to be replaced by sum over Matsubara frequencies. 
Therefore, we do not repeat all the steps presented there, rather, we limit ourselves to convey the relevant results. 

First, we note that in the DE approach applied to this case, keeping up to two derivatives, the Casimir free energy 
can be written as follows: 



(18) 



where the two dimensionless functions bo and hi can be obtained from the knowledge of the Casimir free energy for 
small departures around the = a = constant case. Indeed, setting: 



V'(xn) = a + T)(x\\) 



one expands Tp in (17) in powers of rj, up to the second order. Thus, 



Tp(a, V ) = Tf\a) + rV>(a, V ) + rf\a,r,) + 



(19) 



(20) 



where the index denotes the order in r\. 



For the expansion above, T g is proportional to the area of the mirrors, L . In terms of the Matsubara frequencies: 



2j ?p, n € Z, the explicit form of the zero order term per unit area is as follows: 



iy(a) 
L d-1 



1 

2p 



E 



(27T) 



P|| , 



1 



^-2ay / wJ +p 2 



a a p 



(21) 



where: 



&o(0 



E 



(27T) 



P H 1 



n=—oo 

- « 2 <«-ir(^=i) ^ 

>• 2 ' 71 — — oc 



dpp d - 2 log 



(22) 
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is the dimensionless function which appears in the DE for the zero order term (£ = a/ j3). 

(2) 

Regarding T ^ , which is necessary in order to find 62, the result can be presented in a more compact form in terms 
of its Fourier space version. Defining the spatial Fourier transform of r\ by: 

^ll) = / ^T elk "' X " ^ k ll) ( 23 ) 

we have 

^? = \J^f m W\W\ 2 (24) 

where 

/ (2) K,k„) = -- £ / T^y/rt + p f + UJ ^ 2 + tow + k ") 2 

' m=—oo J y ' 



1 - exp ( - 2aJujl + pj) exp [2a^(w m + oj n ) 2 + (py + ky) 2 ] - 1 



a -<^)jrf»(- ;n)a | k|| |) (25) 



with 



f d a_1 Dn r / / 

F(%;n,|l|||) = -2£ ]T / (2^H V^O 2 + p| + n) ^ + (p " + ^ 

x 1 = p 1 r ]■ , (26) 

1 - exp [ - 2 /(27rm<£) 2 +pf] ex P { 1y/[M™ + n ) C] 2 + (P|| + 1||) 2 } - 1 J 



which is also a dimensionless function with dimensionless arguments. We have made explicit the fact that the result 
will only depend on the modulus of 1|| , as any dependence on its direction may be got rid off by a redefinition of the 
integration variables. 

The coefficient 62 can be immediately defined in terms of F^ 2 K Indeed, 

1 r 0F( 2 >(£;n,|l|||)-i , . 

"<«-*f w UmM" (27) 
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In the following subsection, we consider the low and high temperature limits of the two coefficients b and b 2 ; since 
they determine the form of the DE in the corresponding limits. We note that the relevant scale to compare the 
temperature with is the inverse of the distance between the mirrors. Thus, in terms of the variable £, the relevant 
cases are: £ — > (zero temperature limit) and £ — > 00 (infinite temperature limit). We discuss them below. 



1. The zero and high temperature limits 

The zero temperature limit corresponds to £ — > 0, and it can be implemented by replacing a sum over discrete 
indices by an integral over a continuous index. Defining fco = 27m£, we get an integral over ko; the Jacobian being 
l/(27r£). The results for the two coefficients, in d dimensions (we introduce the number of dimensions as an explicit 
parameter), are: 



N*0] c «i~^/^iofi 



1 - e 



-2|pi|| 



= 60(d) 



(28) 



and 



\b (d (* - 1 \ dF ti d iM 



62(d) 



(29) 
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6 2 (d) 
60(d) 




d= 1 




0.435 


d = 2 


1+6C(3) 
12C(3) 


0.569 


d = 3 


2/3 


0.667 


d = 4 


-C(3) + 10C(5) 
12C(5) 


0.737 


d = 5 


10^-21 

H>7r2 


0.787 


d = 6 


-2C(5) + 7C(7) 
6C(7) 


0.824 



TABLE I. Values of the ratios 



''2 1''! 
60(d) 



for the lowest dimensions. 




FIG. 1. The function 62/60 tends to 1 for large values of d. 



where: 



Ff\d,\l {l \) = -2 J ^IPiillP 



"'ill- 



,-2\ Pll \ e 2| P ||+/||| 



(30) 



It is possible to give a more explicit expression for 6 2 (d) , since this coefficient may be obtained by taking derivatives 
inside of the integrand of (30). Also, the bo(d) coefficient can be exactly evaluated as a function of d. We present, in 



the following table, the ratio between the two coefficients as a function of d, for 1 < d < 6: 

It is interesting to remark that the relative weight of the NTLO correction grows with the number of dimensions. 
Indeed, the general results for ba(d) and 62 (d) in an arbitrary number of dimensions are: 



and 



6a (<0 



bo(d) = - 



1 7r d / 2 



r((d + i)/2) c(d + i) 

(4 7r )(d+l)/2 

(d-3)(d - 1) 



(31) 



12tt 2 2 d 

+ 7r 3 / 2 - d (d + i)r 



d 

d+l 



r ( 2 - - ) c(2 - d) 



C(d+l) 



(32) 



These expressions are consistent with those derived in Ref. [10] . using a different method and in the context of the 
dynamical Casimir effect. Figure [l] shows that the ratio b 2 (d) /bo{d) is an increasing function of d, tending to 1 as 
d — » 00. 

Let us now consider the very high (infinite) temperature limit. When £ >> 1, we see that only the n = term in 
the sum representing 6 yields a non-vanishing contribution, 



(27T) 



d-l 



log 1 



-2|pnh 



(33) 



or, introducing explicitly the dependence on the number of space dimensions, d, 

[&o(&d)]e»i ~ e[&o(£,d-l)k->o = 6o(d-l) 



(34) 



a reflection of the well known 'dimensional reduction' phenomenon at high temperatures, for bosonic degrees of 
freedom. 

For the 62 coefficient, a similar analysis shows that only m = has to be kept, and: 



£[6 2 (£,d-l)] 



b 2 (d-l) 



(35) 



Putting together ( 34 ) and (|35j) , we finally get for the DE up to the second order in the high temperature limit 

(dip) 2 



In particular, the free energy reads, for d = 3, 

[l>ty,3)],,_ ~- I ^ 



1 



[V(X|| 



Id— 1 



6a(d-i; 



[V>(xi|)] 



d-1 } 



\i/>/P»l 



[^(x||)] ; 



{1 + 0.569 (av) 2 } 



(36) 



(37) 



As an example, let us now apply the results above to the evaluation of the Casimir interaction between a sphere 
and a plane at very high temperatures. The sphere has radius R, and is in front of a plane at a minimum distance 
a (a <C R). Although the surface of the sphere cannot be covered by a single function z = ip(x\\), as in previous 
works, we will nevertheless consider just the region of the sphere which is closer to the plane [I]. We shall see that this 
procedure still produces results which are quantitatively adequate within the present approximation and assumptions, 
even beyond the lowest order. 

The function tp is 



R l-i/l- 



R 2 



(38) 



where we are using polar coordinates (p, (p) for the x 3 = plane. This function describes the hemisphere when 
< p < R. The DE will be well defined if we restrict the integrations to the region < p < Pm < R- We will assume 
that pm/R = 0(1) < 1. Inserting this expression for tp into Eq.(37l and performing explicitly the integrations we 
obtain 



[I>(V,3)] 



C(3)fl (-. 
8/3a V 



0.569|log(| 



(39) 



Note that, as long as a <C R, the force will not depend on pm- As expected on dimensional grounds, the R/a 2 
behavior of the leading contribution in the zero temperature case changes to R/a/3 at very high temperatures. The 
NTLO correction is analytic when written in terms of derivatives of the function tp, but non-analytic in j|. This 
behavior has been already noted in numerical estimations of the Casimir interaction between a sphere and a plane in 
the infinite temperature limit, for the electromagnetic case in Ref.|llj. 

It is interesting to remark that the expression for the free energy at high temperatures in d = 3 is quite similar to 
that corresponding to the electrostatic force F z between two surfaces held at a constant potential difference V [7] 



F, 



d 2 x. 



1 



(40) 



Therefore, when considering an arbitrary surface over a plane, the high temperature limit of the free energy will have 
the same behavior than the electrostatic force. For instance, from the results of 0, for a cylinder of radius R and 

length I at a distance a of a plane, we expect the leading term of the free energy to be proportional to h\ ^5 while 



its NTLO correction must be a coefficient times 

Going back to the general case, at intermediate temperatures the coefficients 62(C) given by Eq. 
Eq.(22 1 should interpolate between their zero and high temperatures values. This is shown in Figs. p2 



27} and 6 (£) of 
and [3j where we 

plot the ratio 62/60 a s a function of the dimcnsionless temperature £ for d = 3 and d = 4, respectively. The plot for 
d = 3 interpolates between the value 6 2 /6o = 0.67 for zero temperature, and 0.57 at high temperatures. On the other 
hand, for d = 4, 62/60 interpolates between 0.74 for zero temperature and 0.67 at high temperatures. These limits are 
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FIG. 2. Ratio between coefficient 62(C) given by Eq.(27) and the coefficient 60(C) °f Eq. ( 22 1, as a function of the dimensionless 
temperature £ for d — 3. The plot interpolates between the value 62/60 = 0.67 for zero temperature, and 0.57 at high 
temperatures. 




0.2 0.4 0.6 0.8 



1 1.2 1.4 1.6 i.e 



FIG. 3. Ratio between coefficient 62(C) given by Eq.(27l and the coefficient 60(C) °f Eq.(22l, as a function of the dimensionless 
temperature C f° r d = 4. The plot interpolates between the value 62/60 = 0.74 for zero temperature, and 0.67 at high 
temperatures. 



in agreement with the results in Table 1. The ratio 62/60 gives a quantitative measure of the relevance of the NTLO 
correction to the PFA. Note that, both for d = 3 and d = 4, it converges quickly to the infinite temperature value. 

Finally, it is worth stressing that, as the coefficients 60 and 62 are functions of £ = Tip, the evaluation of these 
functions is crucial in order to compute the Casimir free energy using the DE in any concrete example, at a fixed 
temperature. The previous plots describe the dependence of their ratio with distance, at a fixed temperature. In Fig. 
[4] we plot &o(£) in 3 + 1 dimensions. We can see that, at low temperatures, the curve is very flat. Indeed, it is well 
known that the low temperature corrections to the free energy for parallel plates are proportional to £ 3 for £ -C 1, 
and this behavior is well reproduced in the numerical evaluation. Moreover, the function 60(C) acquires very quickly 
the linear behavior expected at very high temperatures. For the sake of completeness, in Fig. |5j we plot the function 
62(C), which has similar characteristics. 

These results may be useful to understand the nontrivial interplay between geometry and temperature for open 
geometries, like the sphere-plate and the cylinder-plate configurations, described in Ref.|12^. Indeed, it was pointed 
out there that local approximation techniques such as the PFA are generically inapplicable at low temperatures. From 
our results we sec that both functions 60 and 62 approach their high temperature behavior for relatively low values of 
£. Therefore it would not be valid to insert the low-£ expansions of these functions into Eq.(18), and then apply the 
result to open geometries for which the condition £ <C 1 is not satisfied. 



10 




FIG. 4. The function feo(C) in 3 + 1 dimensions. 




FIG. 5. The function 62(C) in 3 + 1 dimensions. 



IV. NEUMANN BOUNDARY CONDITIONS 



Again we start from the functional representation of Zp{ip) given in @, but now we use the functional 8 functions 
which impose Neumann (rather than Dirichlet) boundary conditions on the two mirrors. We assume the mirrors to 
be characterized by the same surfaces we used in the Dirichlet case. 

The boundary conditions may be written as follows: 

L) [d d (p{x\\,x d )] Xd=a = 

R) [d n (p(x\\,x d )] Xd=i ,( Xll ) = , (41) 
where d n — n p <9 M , with the unit normal to the surface x d — VK^il) = ; 

JV"(s„)=#-ag0 o #C|,), (42) 



and \N(x\\)\ = ■ s /g{x\\). 
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To exponentiate the <5-functionals, we again introduce two auxiliary fields, Xl and A^: 

S L (ip) = J T>\ L e i f ddx tt XL( - x "' )[dd ' p( - x n' Xd ' >]:c d=" 

S R ((p) = J V\ R e^ ddx]] V^iO x r( x iO l&n<p(x\\,xd)]x d =i,<.x\\) 

= J V\ R e^ ddx]] Aj?(x " ) { dN ^ x \\< x ^*a=i>(*\\) (43) 
where we introduced the notation <9jv = N^d^. Thus, using those exponential representations, we derive: 

Z p (%l>) = Jv(pVX L VX R e-SoM + ifd^xJNixMx) ) ^ 

where, by analogy with the Dirichlet case, we introduce the 'current' Jn(x): 

Jn{x) = X L (x\\) d d S(x d ) + X R (x\\) d N 5(x d - ip(x\\)) . (45) 

Note that there is no need to get rid now of any metric-dependent factor, as we did for the Dirichlet case. 
The integral over <p becomes then: 

Zp(ip) = Zy VX L VX R e 11 II , (46) 



where: 



f d a h 
Ull(x 11 ,x ]1 ) = - J 



V LR (x h x\, 



ddk w pifeii-^H-xj,) N 

) d 2 
d ^ ,*,-(-,-' 1 ).-|*,|*(-' 1 )NI +< *l|-^ !e fi) 



(27r) d 

f dd h 

U R l(x 11 ,x ]] ) = - J j—y d 



e ife r (x||-X||) e -\l 



2 



\} RR {x h x\^ = f e **ii-(«ii-»M) e-lfciiK^ii)-^^)] 



X 2 



The free energy Tpty) is then 



(27r) d 

\{ - |*„ | -ik r [d^{x\\)+ d^{x[)] + ±-(k r d^{x\\)k r d^{x[))) (47) 



iyvo - ^ Trl °s u ' ( 48 ) 



which, as in the Dirichlet case, does contain 'self-energy' contributions, to be discarded here by the same reason as 
there. Again, <]/ is assumed to be time independent. 

A. Derivative expansion 

Assuming that one could proceed as in the Dirichlet case, keeping up to two derivatives, the derivative expanded 
Casimir free energy could be written as follows: 

with two new dimensionlcss functions Co and C2. Those coefficients may be determined from the knowledge of the 
Neumann Casimir free energy for small departures around the V(x||) = a = constant case, up to the second order in 
the departure. Again, we focus on the cases of purely quantum or purely thermal effects, except for the more realistic 
case of d = 3. As we will show in what follows, the NTLO term is quadratic, except when d = 2 at zero temperature, 
or when d = 3 and there is a non zero (finite or infinite) temperature. 
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We present the calculation of the terms contributing to that expansion in Appendix A. 

It is quite straightforward to see that the zero order term coincides with the one for the Dirichlet case; namely: 
c = b . 



(2) 

rjj ; has the following form: 



where 



with 



,(2) _ 



+ 00 



1 f d^k. 



2 J (2tt) c 



L fl ( 2 )(0,k„)|^(k,|)| 



5 ( 2 )( Wn> k|,) = -l £ 



d d 1 p\\ [w m (w m +UJ n ) +P|| • (P|| +k||)] 2 

^) d ~ l . /^ + p| v /K l +^) 2 + (P||+k||) 2 
1 



1 - cxp ( - 2ay / w| + pjj) ex P [2a-v/( w ™ + w «) 2 + (P|| + k ||) 2 ] - 1 
= a - (d+ 2) G ( 2 )(« ;n)a | k||D 



+ ■00 



G( 2 )(C;n,|l|||) = -2e £ 



d Pll 



[(2^) 2 m(m + n) + p,| ' (Pll +1||)] 2 



1 



(2ntfm* + p 2 v /(2 7 r0 2 (m + n)2 + (p ||+ l || 

1 



1 - exp ( - 2^/(27r£) 2 m 2 + p 2 ) exp [2y/{2n^(m + n) 2 + (p,, - 1 

P2(0 



2 L a|i|i| 2 



n-vO,|l||KO 



(50) 



(51) 



(52) 



(53) 



1. The zero and high temperature limits 

As before, the zero temperature limit can be implemented by replacing a sum over discrete indices by an integral 
over a continuous index. The result for the coefficient ci in d dimensions is: 

[ c a( d '0] € «i = [ ^|^ia l ^ l]'ii-»° = c a( < » ( 54 ) 

where: 



4 2) (^ud=-2 /|ah-(fn+y] * * . (55) 

v ' l||ly y (2?r) d |p||| + 2||| 1 - e - 2 IPiil e 2 l p n +z iil - 1 

For d — 1, the coefficient C2 coincides with its Dirichlet counterpart 62- In higher dimensions, the structure of the 
form factor is different. We present, in Table 2, the ratio between 02(d) and co(d) = &o(<i) as a function of d, for d ^ 2: 

We see that the ratio c 2 (d)/c (e?) is non-monotonous and always negative for d 7^ 1. We have also checked that 
C2 (d)/co (d) — 1 for large values of d. 

In the particular case d — 2 is not possible to compute the coefficient by introducing the derivative with respect to 

| 2 inside the integral in Eq.(|55 [), because of infrared divergences. This is a signal of a branch cut at zero momentum, 
as we will show in Section II V Bl 

The high temperature limit can be obtained, as for the Dirichlet case, taking the limit £ >> 1. "Dimensional 
reduction" takes place, and the free energy becomes 

MiMWm ~jif dd S\ {bo(d - 1) ^ja=i + ^ i) [^ff -1 } ■ ( 56 ) 

Regarding intermediate temperatures, Fig. [6] shows the ratio between coefficients c 2 (£) and co(£), as a function of 
the dimensionless temperature £, for the Neumann boundary condition in d = 4 dimensions. The plot interpolates 
between the value C2/C0 = —1.00 and -1.36 at zero and high temperatures respectively, in agreement with the results 
shown in Table 2. 
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c 2 jd) 
cold) 




d= 1 


i ^ t n 'i 
3 3C(2) 


0.435 


d = 3 


3 3C(4) 


-1.360 


d = 4 


6 12C(5) 


-1.002 


d = 5 


1 9C(4) 
5C(6) 


-0.915 


d = 6 


7 2C(5) 
6 C(7) 


-0.890 


d = 7 


4 4(,0'M 

3 C(8) 


-0.886 



TABLE II. Values of the ratios 



c 2 (d) 
co(d) 



for the lowest dimensions. Note that we have excluded the case d — 2. 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 



FIG. 6. Ratio between coefficient C2(£) and the coefficient co(f;), as a function of the dimensionless temperature £, for the 
Neumann boundary condition in d = 4 dimensions. The plot interpolates between the value C2/C0 = —1.00 and —1.35 at zero 
and high temperatures respectively. 

B. Non analytic terms: d = 2 with T = 0, and d = 3 with T > 



Let us consider the particular case of d = 2 at zero temperature. As shown in Appendix A, for small departures of 
the plane-plane geometry V'( x ||) = a + 7 7( x ||)) the correction to the Casimir energy reads, up to second order in r) 



r (2) = 



1 



d 2 fc|| 

(2^ 



[» w (*i)]^l«Wf. 



with: 



5 (2) W = -2 



rffan h • fan +%)]' 



(57) 



(58) 



(2tt) 2 byl by + fc||| l - e _2a l p iil e 2a|p ii +fc » l - 1 ' 

Naively, one would expect the form factor </ 2 )(fcy) to admit an expansion in powers of fc 2 , which is the necessary 
condition in Fourier space to produce a DE in configuration space. However, this is not the case for d = 2, as suggested 
by the fact that the formal expression 



dg 



(2) 



-P\\ 



1 



d 



<VA-if lfc "~ f ° 2 ./ (2tt) 2 |p„| (1 - e -2a| P|l |) dkf. 



1 



[p\\ ■ fall + fc ll)] 2 _ 

fall I fay + fcyl e 2 "IPii+ fe iil - 1 



(59) 



J fen — ^0 



has an infrared logarithmic divergence at pi =0. 

The behavior of </ 2 )(fc||) for small values of k\\ can be determined by studying the integral that defines it in Eq.(58) 
in the region a\p» \ <C 1. In this region, and assuming also that a\k» \ <C 1 one can make the approximation 
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and compute the integrals analytically. In this way, it is possible to show that 



/ 2 )(fc,|) « g^(0) - log^jja 2 ) + 0(k 2 /a 2 ) . (61) 



This behavior of g^ 2 \ that we confirmed with a numerical evaluation of Eq.(58), shows that a local DE breaks down for 
Neumann boundary conditions at d = 2. However, one can still perform an expansion for smooth surfaces, including 
nonlocal contributions in the Casimir energy. For instance, in the present case, the NTLO correction to the PEA will 
be nonlocal and proportional to 

d 2 X||?7(a;||)V|log(-a 2 V 2 )77(a;||) . (62) 

As we will describe more generally in the next Section, the breakdown of the local expansions is related to the 
existence of massless modes in the theory. These modes are generally allowed by Neumann but not for Dirichlct 
boundary conditions, that impose a mass gap of order 1/a. 

The logarithmic behavior of the form factor in d = 2 induces a similar non-analiticity for d = 3 at finite temperature. 



Indeed, the m — n = term in the finite temperature form factor given in Eq.(52) is formally identical to the Neumann 
form factor in d = 2. Therefore, in an expansion for small values of \k\\\, in addition to a term proportional to 
A: 2 , there is a contribution proportional to (Ta)k 2 log(/c 2 a 2 ) at any non-vanishing temperature, which is not cancelled 
by the rest of the sum over Matsubara frequencies. 

V. HIGHER ORDER TERMS IN THE DE 

In this Section we discuss some general aspects of the derivative expansion, related with the calculation of higher 
orders and the eventual breakdown of the expansion. 

In this and previous works we considered just the NTLO correction to the PFA, which contains up to two derivatives 
of the function -0 in the free energy. We expect the next to NTLO (NNTLO) order corrections to contain terms of 
the form 

|<9V| 4 , VWI 2 ^, ^ 2 d 2 ^d 2 ^, i/} 2 d a d^d a d ip, ^d 2 d 2 ^ , (63) 

and terms containing more derivatives for higher orders. 

The main question to be answered is whether the free energy can be expanded or not in local terms up to any 
desired order. For this to hold true, a necessary condition is that the expansion must hold true for a particular case: 
when ip = a + rj, with rj -C 1, and one keeps just the quadratic term in rj. 

For simplicity, we deal with the T — case in d = 3 dimensions, an example that will illuminate several aspects of 
the problem. From the previous sections, we see that the quadratic contribution to the Casimir energy has the form: 

r£) = i?J ^^ (2) (°' ak ii)l^ k ii)l 2 ' ( 64 ) 

where the form factor depends on the boundary conditions of the quantum field. In most of the paper we considered 
the case of static surfaces, but here it will be useful to discuss the more general case in which the right mirror can be 
in motion. In this situation, on general grounds, we expect the form factor to be a function of a{k^ + k 2 ) 1 / 2 = a\k\\ |, 
and of course the explicit calculations confirm this fact. Within these approximations, the Casimir energy will not 
admit an expansion in derivatives if the form factor includes, for instance, odd powers or logs of its argument. Note 
that at higher orders, and relaxing the condition of a quadratic approximation in rj new non-analytic terms may arise, 
which can be of the same order in the DE as the ones that come from the term quadratic in rj. 

In 3+1 dimensions, this question can be answered from the explicit expressions of the form factors |13| presented in 
Appendix B. For Dirichlet boundary conditions, the expansion of the form factor contains, in addition to even powers 



of the argument, a term proportional to a |fcii| (see Eq.(82)). The non-analytic term becomes a nonlocal contribution 



in configuration space, that cannot be expanded in derivatives of r\. Note that this contribution does not depend on 
the distance between mirrors. There is a simple interpretation of this term: when considering a flat moving boundary, 
photon creation produce an imaginary part in the vacuum persistence amplitude when rotated from Euclidean to 
Minkowski spacetime. Therefore, the expression for T^} cannot be analytic in fcjj. For a single nonrelativistic mirror, 
this will lead to a dissipative force proportional to the fifth time-derivative of the position. Indeed, the proper analytic 
continuation of the Euclidean term proportional to |fco| 5 to real time, produces this dissipative force. For the Dirichlet 
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case, there are no non-analyticities dependent on the distance between mirrors. Physically, this is due to the fact 
that there is a frequency threshold to produce real photons between mirrors, which is of order 1/a. The conclusion 
is that, for Dirichlet boundary conditions, the a— dependent part of the DE is well defined up to any order. On the 
other hand, the a-independent part contains a non-analytic contribution related to the possibility of creating photons 
of arbitrary low energies from the vacuum. 

The situation changes for the case of Neumann boundary condition. As shown in Eq. ( 84 ) , in addition to the non- 
analytic term proportional to |A;||| 5 a 5 , which as before produce a contribution to T^} that is independent of a, there 
is a term proportional to \kn | 3 a 3 . The physical reason of the existence of this term is again clear using the connection 
with the DCE. Indeed, this term comes from the possibility of creating TM photons between mirrors, moving parallel 
to the mirrors, for which there is no energy threshold. If the DE is used to compute the force on a moving mirror, 
this a-dependent dissipative contribution to the force will be missed, i.e. it is only possible to get the dispersive part 
of the force. On the other hand, for a static and non-fiat mirror, this term produces a nonlocal component of the 
force, which will be smaller than the NTLO, but the dominant part of the NNTLO. 

It is interesting to remark that both the non-analytic contributions proportional to \k\\ | 5 for Dirichlet and Neumann 
boundary conditions in 3 + 1 dimensions, and the |fc||| 3 term for the Neumann case, could be derived from the form 
factors described in the previous section by analyzing the infrared behavior of the integrals in p» , as we did in Sec. 



IV B Moreover, it is clear that all of them have the same physical origin: the existence of massless degrees of freedom. 



For Neumann boundary conditions at T = 0, the non-analyticities show up in the NTLO for d = 2, and in the NNTLO 
for d = 3. For Dirichlet boundary conditions the non-analytic term is independent of the distance between mirrors, 
and therefore does not contribute to the Casimir force. 

The situation is analogous to that of effective field theories that involve massless particles. In that case, in addition 
to local terms in the effective action, there are nonlocal (or non-analytic contributions) which can be interpreted as 
arising from the fact that there is no threshold for creating such particles [5] . A prototypical example is quantum field 
theory under the influence of external (classical) backgrounds. For massive quantum fields in curved spaces |14j . the 
effective action and the energy momentum tensor of the quantum fields can be approximated by a DE (usually known 
as the Schwinger DeWitt expansion in that context). Each subsequent term in the expansion contains additional 
derivatives of the metric and inverse powers of the mass of the quantum field. The expansion is valid as long as the 
typical scale of variation of the classical background is much larger than the inverse mass. However, for massless 
quantum fields, it is necessary to consider nonlocal contributions. In our case, the role of the background is played 
by the curved surface, and there are both massive and massless excitations: the massive ones are the Dirichlet modes 
(TE photons) inside the " cavity" . The massless ones are the Neumann modes (TM photons) with momentum in the 
direction parallel to the plates, and TE and TM photons outside the cavity. 

The physical picture suggests possible ways out to improve the PEA even beyond the NTLO correction. This would 
involve a separate treatment of massless and massive excitations. We hope to address this issue in a future work. 



VI. CONCLUSIONS 



We have obtained explicit expressions for the NTLO term in a DE for the Casimir free energy for a real scalar field in 
d spatial dimensions. The field satisfies either Dirichlet or Neumann boundary conditions on two static mirrors, one of 
them flat and located at the Xd — plane, while the other is described by the equation: — ip{x\, X2, ■ ■ ■ , Xd~i)- We 
have shown that, for Dirichlet boundary conditions, the NTLO term in the Casimir energy is always of quadratic order 
in derivatives, regardless of the number of dimensions. Therefore it is local, and determined by a single coefficient. 
We evaluated the ratio between that coefficient and the one for the PEA term, for different values of d at zero and 
high temperatures. 

We have also shown that the same holds true, if d ^ 2, for a field which satisfies Neumann conditions. When 
d = 2, the NTLO term becomes nonlocal in coordinate space, which is a clear manifestation of the existence of gapless 
excitations allowed by the Neumann conditions |15j . It may be seen that among all the possible combination of linear 
boundary conditions on the mirrors, just this case, Neumann conditions on both mirrors, can produce these modes. 

When including thermal effects, we have shown that, for Dirichlet mirrors, the NTLO term in the free energy is also 
well defined (local) for any temperature T. Besides, it interpolates between the proper limits: namely, when T — > 
it tends to the one we had calculated for the Casimir energy, while for T — > oo it corresponds to the one for a d = 2 
theory, realizing the expected dimensional reduction at high temperatures. On the contrary, for Neumann mirrors in 
d = 3, we found a nonlocal NTLO term for any T > 0, which vanishes linearly when T — > 0. This leaves room, when 
the temperature is sufficiently low, to use just the local term (of second order in derivatives) as the main correction to 
the PFA. But of course, the nonlocal term will always break down for a higher temperature, whose value will depend 
on the actual shape of the surface involved. We stress once more that this non-analytic behavior is a consequence of 
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the Neumann boundary conditions, and may not be present for imperfect boundary conditions, as those considered 
in Ref.0. 

In the course of our derivations we have obtained integral expressions for the momentum space kernels which 
determine the quadratic contribution to the free energy for small departure from the planar case. Those kernels are 
well defined in any number of spatial dimensions and temperatures, and agree with the known results for d = 3 
and T = [13] . They can be used to extract the NTLO terms, be they local or nonlocal. Although for the static 
cases we have considered in this article they are only needed for time independent configurations, we also present the 
expressions for the kernels at non-zero frequencies. 

APPENDIX A: EXPANSION TO ORDER 0(?7 2 ), NEUMANN CASE. 

We present here the main steps and intermediate results corresponding to the calculation of the free energy up to the 
second order in the function 77, which measures the departure from the planar case. We assume that ip(x\\ ) = a + rj(x\\), 
with a equaling the average of tp. Namely, we want to construct the terms in 

r> = rj°> + r« +rf ) + ... (65) 

where the index denotes the order in rj. The term of order 1 vanishes, and, in terms of the expanded matrix elements 
of U, we may write the more explicit expressions: 



where: 





p(0) 




r (2) 


(2,1) 
~~ 


— Tr[( 

2(3 V 


(2,2) 

~ 





2/3 

r (2,l) + p (2,2) ^ (66) 



( U (0))- 1 U (2)" 

(U^rV^U^r'uWj . (67) 



In order to simplify, and at the same time render the expressions more compact, we shall keep the component of 
momenta to appear below continuous, as if they corresponded to zero temperature. In order to obtain the proper 
finite temperature expressions one should just replace integrals over the component of the momenta by sums over 
Matsubara frequencies. Namely, 

f ^... A(fc0j . g A(w n , ...) . (68) 

n— — oo 

Let us consider the explicit form of (U^ -*) and its inverse, since both of them are required to construct the terms 
contributing to Tp above. We first note that the zero order term is given by: 

TJ (o) r M _ ( d d d' d A{x - x')\ Xd . x ^ d d d' d A{x-x')\ Xd ^, x , d ^ a \ 

V ^ X \\> - {d d d' d A(x-x')\ Xd ^ x ^ d d d' d A(x - x')\ XM ^ a J ' ^> 



which, because of its independence of 77, can be conveniently Fourier transformed in the parallel coordinates: 

d d k\\ 
(2?rp 



U(°)(x|„xf,) = U<°>(s„ -x' {l ) = [ *^ e <VC*ii-*f|) U(°)(fe„) , (70) 



where 



Thus: 
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Regarding U' 1 ), we see that 



< = 4 X 2 = o 



(73) 



and 



i /• d d /c 



2 7 (2tt) 



Jl.e w n-(«ii-*n)[|jfe||| a » 7 (a;||)+tJfc a ft,tj(x||)]e-l fc iil 



(74) 



(2) 

Finally, to the second order, = 0, and 



u s?>ih4) 



1 f d d kn 



Jk r (x«-x' u ) 



k a k 



2 J {2ir) d 



(|fc|||3 + W dad ' fj) v{X{])v(x '" ] 



(75) 



(2) (2) 

Regarding and U^jj, they are non- vanishing, but it may be seen that they do not contribute to the second order 
term. 
Thus, 



^(2,1) _ 



2/3 



Tr 



= j ddx \\ ddx '\\ {<1)~\h^rr^\\) 

and (using the properties of the matrix elements under the exchange of arguments) : 



(76) 



.(2,2) 



2/3 



Tr 



Kir^(<ru (i) 



RL 



2/3 



Tr 



LR 



(77) 



In Fourier space, after some algebra, one then finds: 



.(2) 



1 [ d d ku 



(2tt) 



(2) (fc||)|^|| 



(78) 



with: 



/ 2) (fc||) = -2 



dd P\\ [P\\ ' (P|| + fc ll)] 1 1 

(2ir) d \p» | \p\ I + fcii | 1 - e^ 2 *" 1 e 2a|p ii +fc ii | - 1 



(79) 



APPENDIX B: EXACT EXPRESSIONS FOR THE FORM FACTORS IN 3 + 1 DIMENSIONS. 

In this Appendix we present exact expressions and series expansions for the form factors /^(^||) an( i fl^ 2 H^||) a ^ 
zero temperature and d — 3. These expressions have been previously obtained in Ref.[13] (see also Ref.[16]). We will 
use the dimensionless quantity x — a|fc|||- 

For Dirichlet boundary conditions the form factor f^> reads 

f < 2 > ffc„ ) - -2 f^L \P\\WP\\+ k \\\ m 

J [ V J (2tt) 3 (1 _ e -2«IPiil)( e 2«IPii+*iil - 1) ' 1 ' 

An explicit evaluation gives [13] : 

x 3 Li 2 (er 2x ) x 2 U 3 (e- 2x ) xLi 4 ( e - 2x ) Li 5 (e~ 2:r ) 7r 2 Li 2 (l - e~ 2x ) 



a 5 f^(x) 



48tt 2 24tt 2 16tt 2 16tt 2 240x 

-U e (e- 2x ) xHog(l~e- 2x ) ir 2 x 

32n 2 x 120tt 2 240' 1 ' 
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where Li n denote Polylogarithm functions. This result can be expanded in powers of x as follows: 

5 (2) _^_tt 2 ^ 2 (45 + TT 4 ) x 4 x 5 (315 - 2ir 4 ) x 6 

a J W » 12Q 108Q + 27000tt 2 720tt 2 + 793800tt 2 



(7T 4 - 105) X 8 (l65 16335 ) X 

5103000tt 2 30240tt 2 



10 

12 



O (x 12 ) (82) 



which shows the presence of a non-analytic term proportional to x , whic h is independent of a 



For Neumann boundary conditions the form factor g^> is given in Eq.(79) with d — 3. Evaluating explicitly this 
integral it is possible to show that [13] 

7 5„(2)^_ 1 , ^„ T5 ./„-^ , f 1 ^ W,„-^ 5xLi 4 (e- 2 -) 7Li 5 (e- 2 -) 



a b 9^{x) = _ _ + 1 xLi 2 e- 2 ^ + - - — Li 3 e 



24V27T 2 J v ; V 16 487r2 / ' 3271 " 2 3271-2 



TT X 



(83) 



* a Li 2 (1 - e- 2 ") | ^ (e- 2 -) ._, 

240a; 32tt 2 x ^ 120tt 2 720 

Expanding this result for i « 1 we obtain 

«, ,,,, x TT 2 (30-tt 2 )x 2 x 3 (1095 + 50tt 2 + tt 4 ) x 4 llx 5 (2205 - 42tt 2 - 2tt 4 ) x 6 

a CT MX) ~ h ^ - h - - h - - h 

y K ' 120 1080 64 27000tt 2 720tt 2 793800tt 2 

(l007r 2 + 7^ 4 - 3045)x 8 / 47 1 ^ Tioks + Is^Wo \ _io . Q( ^ (M) 

35721000tt 2 \ y 381024007T 2 23328000 246985200 128tt 2 J 

We see that the Neumann form factor has non-analytic terms proportional to x 3 and x 5 . 
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